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Abstract: We consider = 4 super Yang-Mills theory on a four-dimensional lattice. 
The lattice formulation under consideration retains one exact supersymmetry at non-zero 
lattice spacing. We show that this feature combined with gauge invariance and the large 
point group symmetry of the lattice theory ensures that the only counterterms that appear 
at any order in perturbation theory correspond to renormalizations of existing terms in 
the bare lattice action. In particular we find that no mass terms are generated at any 
finite order of perturbation theory. We calculate these renormalizations by examining the 
fermion and auxiliary boson self energies at one loop and find that they all exhibit a common 
logarithmic divergence which can be absorbed by a single wavefunction renormalization. 
This finding implies that at one loop only a fine tuning of the finite parts is required to 
regain full supersymmetry in the continuum limit. 
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1. Introduction 

A/" = 4 super symmetric Yang-Mills (SYM) theory in four dimensions is both a fascinating 
quantum field theory in its own right and in addition plays a crucial role in the well known 
AdS/CFT correspondence ||l|, ^ where it is thought to provide a dual description of type 
IIB string theory in AdSrj x 5*5 space. 
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The strong coupling, large- limit {N being the number of colors) of this gauge the- 
ory has been extensively studied since the corresponding dual theory reduces to a weakly 
coupled supergravity theory, which describes the low energy limit of type IIB string theory. 
However it has proven difficult to go beyond this leading supergravity approximation in 
most situations. Having a lattice formulation of the super Yang-Mills theory would clearly 
be very advantageous as it would give a non-perturbative definition of the gauge theory 
and offer up new tools to investigate its strong coupling dynamics and thence the dynamics 
of the dual theory. Indeed such a lattice construction would allow for a systematic study 
of the classical and quantum string corrections to the supergravity solution. Unfortunately 
straightforward attempts to discretize the continuum theory are well known to break su- 
persymmetry completely leading to a profusion of supersymmetry violating counterterms 
(four or six counterterms, depending on the gauge group, in the effective action whose 
couplings must be fine tuned to approach the correct continuum limit - see Ref. Q). 

However, recent formulations of super symmetric lattice theories, which retain exact 
supersymmetry at non zero lattice spacing offer the hope of at least partially evading these 
fine tuning problems - the existence of an exact lattice supersymmetry protecting the theory 
from many of these dangerous counterterms 0, |5|, ^, ^, IC , 11, |l^, [l^, 14 ] . See also 
the recent reviews |l6| for further references. In addition, there has recently been a 
great deal of other work, both theoretical and numerical, focused on discrete formulations 
of A/" = 4 SYM see Refs. 0, |l8|, ||, H, |l|, H, ||]. These alternative approaches should 
be viewed as complementary to the lattice construction described in this paper. 

In the case of A/" = 4 SYM the corresponding lattice theory retains only one out of the 
sixteen continuum super symmetric invariances and the question of how much fine tuning is 
required to take a continuum limit of this lattice theory targeting the usual = 4 theory 
has been unclear up to this point in time. This paper aims to address this issue using both 
general arguments valid to all orders in perturbation theory and an explicit calculation of 
the renormalization of the lattice theory to one loop order. 

We will argue quite generally that the symmetries of the lattice theory strongly con- 
strain the possible counterterms that can arise as a result of quantum corrections; we find 
that the only relevant operators that can be induced via radiative effects correspond to 
renormalizations of four marginal operators already present in the tree level theory. These 
operators correspond to kinetic terms and we show, using a topological argument based on 
the exact (twisted) lattice supersymmetry, that no mass terms are induced to all orders in 
perturbation theory. 

The remaining fine tuning question then hinges on what divergences can arise in the 
renormalization of these four bare couplings. We proceed to calculate these divergences 
at one loop using lattice perturbation theory. Exact lattice supersymmetry allows us to 
extract these leading divergences by examining the renormalization of the three types of 
twisted fermion propagator and a single propagator for an auxiliary bosonic field. We show 
that all these exhibit a common logarithmic divergence at one loop. The appearance of a 
single logarithmic divergence ensures that at one loop only finite parts need to be fine tuned 
in order to regain full supersymmetry in the continuum limit. This is a huge advantage of 
this approach as compared to earlier efforts at constructing supersymmetric lattice theories 
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in four dimensions. 

We start with a discussion of the approach to super symmetric lattices through dis- 
cretization of a topologicahy twisted formulation of the super Yang-Mills theory, write 
down the lattice theory and discuss the constraints on its renormalization implied by the 
lattice symmetries. We then derive the Feynman rules governing the perturbative structure 
of the lattice theory and write down the diagrams needed to renormalize the theory at one 
loop. We compute the partition function at one loop and show that it is independent of 
any background fields and furthermore that this is true to all orders in perturbation theory. 
We then evaluate the one loop diagrams, extract their leading logarithmic divergences and 
compute the required renormalization of the lattice theory. We conclude with a summary 
of our main results. 

2. Twisted J\f = A, d = 4 Super Yang Mills 

The key idea that allows us to construct a supersymmetric lattice theory that targets A/" = 4 
super Yang-Mills in the (naive) continuum limit is called topological twisting - see ||l^ and 
references therein^. This twisting process, which in flat space can be thought of as merely 
a change of variables, exposes a scalar, nilpotent supercharge Q. It is the supersymmetry 
associated with this supercharge which can be implemented exactly in the lattice theory. 

To understand how the twisted theory is constructed one needs to examine the relevant 
global symmetries of the continuum Yang- Mills theory. The Euclidean version of the = 4 
SYM theory on can be obtained by dimensionally reducing J\f = 1 SYM theory on M^'^ 
down to M^. The ten-dimensional theory possesses an SO{W) Euclidean (Lorentz) rotation 
group. After dimensional reduction it reduces to 

SO{10)e SO{4)e X SO{6)i , 

where SO{4)e ~ SU{2) x SU{2) is the four-dimensional Euclidean (Lorentz) symmetry on 
and 50(6)/ ~ SU{4:)r is the global internal i?-symmetry group of the dimensionally 
reduced theory. 

The basic idea of twisting is to decompose the fields of the theory in irreducible rep- 
resentations of a twisted rotation group which involves both the usual rotations and the 
R-symmetry^. The global /^-symmetry group of the dimensionally reduced theory contains 
a subgroup SO{4:)ji x U{1). To construct the relevant twist needed for the lattice theory, 
we adopt as twisted rotation group 50(4)' the diagonal subgroup of S0{4:)e x 50 (4) /j. 
The global symmetry is now given by 

G' = SO{Ay X [/(I) 

~ SU{2y X SU{2y X [7(1) 

^This approach to supersymmetric lattices has been shown to be entirely equivalent to the orb- 
ifold/deconstruction formulations pioneered by Kaplan, Unsal and collaborators jl^] which predated the 
twisted constructions in the case of gauge theories. However for simplicity we will mostly be using the 
language of twisting in this paper. 

^The A/" = 4 theory may be twisted in three inequivalent ways ps] , pif however only one of these - 
that due to Marcus and described in the text is ultimately compatible with discretization. 
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Q(3) 
^(3) 



Q(4) 
^(4) 



Notice that since the U{1) part of the global internal symmetry group is undisturbed, it 
remains as a global i2-symmetry of the twisted theory. The supercharges and fermions 
transform under the new rotation group as 

5[/(2)' X SU{2y X Uil) ^ (1, 1) 1 © (2, 2)_i © [(3, 1) © (1, 3)] i © (2, 2)_i © (1, 1) i , (2.1) 

or equivalently 

50(4)' X ^7(1) ^ li © 4_i © 6i © 4_i © li . (2.2) 

2 2 2 2 2 

As a result of this choice of embedding, the twisted theory contains supersymmetries and 
fermions in integer spin representations. They transform as scalars, vectors and higher 
rank p-form tensors: 

Supercharges: Q(°) © Q(^) © Q^^) 
Fermions: © ^^^^ © ^^'^'^ 

We parametrize the fermionic content of the theory by 

7] 1 

The four gauge bosons transform as (2, 2) under the twisted rotation group. We label them 
as a vector field A^. Similarly, four of the six scalars of the theory are now elevated to the 
same footing as the gauge bosons; they also transform as (2, 2) under the twisted rotation 
group. We label them as a vector field B^. The two other scalars remain as singlets under 
the twisted rotation group. We label them by 
theory transform as: 

su{2y X su{2y x u{i) ^ (i, i)i © (2, 2)0 © (2, 2)0 © (1, 1) 

or equivalently 

50(4)' X U{1) ^ li © 4o © 4o - 
We parametrize the bosonic content of the theory by 

1 

4 



4 
6 

4 
1 



(2.3) 



and (h. Thus the bosons of the twisted 



-1 I 



1_ 



(2.4) 
(2.5) 



Ap 



(2.6) 



For future reference the fields of the twisted theory have the following mass dimensions: 



Twisted field 


Ap 


Bp 


</> 




Xpu 


i^P 


i^P 


1] 


77 


Mass dimension 


1 


1 


1 


1 


3/2 


3/2 


3/2 


3/2 


3/2 



In this paper we are interested in Q^^\ the scalar supercharge which is nilpotent 
(q(o))2 . = 0. It is also interesting to see that with this particular choice of twist, the 
resulting twisted fermions are just sufficient to saturate a single Dirac-Kahler field 0. 
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3. Continuum Action and Nilpotent Super symmetry 



Since the two vector fields and of tfie twisted theory transform the same way under 
the twisted rotation group we can describe the theory in a compact way if we combine the 
vector fields into a complex vector field Aa ||26[|^: 



A^ = A^ + iBf,, (3.1) 
Af, = A^- iB^ . (3.2) 

Using these connections one can now define three covariant derivatives and field strengths: 

i^M • = d^, + ig[A^, ■ ], F^, = -'-[D^.D,] , (3.3) 

Vf,- = df, + ig[A^,, ■ ], = -^[V^,V^] , (3.4) 

V^- = d^ + ig[Af„ ■ ], T^^ = -^V^,V^] . (3.5) 

To make contact with the final lattice construction it is useful to assemble the complexified 
gauge fields and the two scalar fields into a single five-component complexified connection 

Aa= (^Af, = Af, + iBf„ A^^s + i^s), a= ,5 ;/i = ,4 (3.6) 

where the fifth component ^5 = (j) and ^5 = (f). Correspondingly the fermions can be 
packaged as five-dimensional scalar, vector and antisymmetric tensors {Tj^^a^Xab)- The 
original twisted theory will then be obtained by simple dimensional reduction of a theory 
in five dimensions. A similar language arises in the orbifold construction of this theory 
where the fermions and bosons transform in representations of SU{5). In addition to these 
fields we introduce one auxiliary bosonic scalar field d for off-shell completion of the scalar 
supersymmetry. 

The nilpotent scalar supersymmetry Q (from now on we denote the scalar supersym- 
metry Q^^^ by Q) now acts on these fields in a simple manner 

QAa = i'a (3.7) 
Q^a = (3.8) 

QAa = (3.9) 

QXab = Tah (3.10) 

Qt] = d (3.11) 
= (3.12) 

The action of the twisted theory can now be expressed in a compact five-dimensional form, 
as a linear combination of Q-exact and Q-closed terms 

5 = QA + SQ.closed , (3.13) 



^Throughout this paper we will be employing an hermitian basis for the generators satisfying Tr (rT*") — 
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where 



and 



A = y Tr (xabJ'ah - -r][Va,Va] + ^vd) , (3.14) 

The invariance of the Q-closed term is a result of the Bianchi identity (or Jacobi identity 
for covariant derivatives) 

1 

eabcdel^cJ^de = Cafecdei^c, Pd, 2?e]] = . (3.16) 

9 

Carrying out the Q-variation and subsequently eliminating the auxiliary field d using the 
equation of motion we can write down the action in terms of the propagating fields. 

\^abcdeXdeT^cXab]- (3-17) 



2 

After a redefinition of the fields grj — )• ry, gijja — )• V'a, 9Xab — ^ Xab and gAa — )• Aa we have 

5 = ^ y Tr ( - \Va,Vi,][Va,Vb] + ]^^a,Vaf " Xab^^a^^b] " rfDa^a 

- T,^abcdeXdeT^cXab]- (3.18) 



2 

The target twisted theory in four dimensions can be obtained by dimensional reduction of 
this theory along the 5th direction. We write down the decomposition of five-dimensional 
fields into four-dimensional fields as follows 

Aa^ A^®(t) (3.19) 

Fab ^ e V^(t) (3.20) 

pa,Va] ^ ^^.,v^] e -\4>A] (3.21) 

^a^^^,®l] (3.22) 

Xab X/iiy ffi Ipf, (3.23) 

The action ( p.l8| ), after dimensional reduction, yields 
5 = ^ y Tr ( - [V^,V,][V^,V,] + ^[V^,V^f + i[^,<^]2 - {v;4>){V^(t>) - x^.^i^.V'.] 

-V'/.^M^ - «V'[<^, ^/.] - - - X^T^fii^u - lpclu\^^ XtJiy]) , (3.24) 

where the last two terms arise from the dimensional reduction of the Q-closed term with 
X*, the Hodge dual of x, defined as x*^u = \^tivp\Xp\ and = \x5p- 

This action can be identified with the twisted A/" = 4 SYM action in four dimensions 



written down by Marcus [26|, up to a trivial rescaling of the fields (with a gauge parameter 
a = 1 in @|). It is important to note that in fiat space, this twisted action is just a 
rewriting of the usual N = A SYM theory in four dimensions and is physically equivalent 
to it. 
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4. Lattice Theory 



Discretization of the twisted theory described in the previous section proceeds straight- 
forwardly; complex continuum gauge fields are represented as complexified Wilson gauge 
links Ua{n) living on links ^u^, a = 1 . . . 5 of a four-dimensional lattice. Since there are 
five such vectors it should be clear that this lattice will have five basis vectors. To ensure 
that the lattice theory enjoys a maximal symmetry we would like these basis vectors to all 
be equivalent. This requirement means that the lattice must possess an point group 
symmetry (the Weyl group of SU{5)). The unique solution to these constraints in four 
dimensions is the so-called lattice. This will hence be the underlying lattice used in our 
work. 

A specific basis for the lattice is given in the form of five lattice vectors 

^ / 1 1 1 1 \ 



/12 V20> 

e, = (0,0,0,--|=). (4.5) 



20> 

These lattice vectors connect the center of a 4-simplex to its five corners. They are related 
to the SU (5) weights of the 5 representation. The unit cell of the A^ lattice is a compound 
of two 4-simplices corresponding to the 5 (formed by the basis vectors e^) and 5 (formed 
by the basis vectors — e^) representations of SU{5). The basis vectors satisfy the relations 

5 1 ^ 

^ ] ~ 0) ■ ~ (^^mn j ^ ^ (^m)/i(^m);/ ~ '^/ii^j A*; ^ ~ 1) ' ' ' > 4. (4.6) 

m=l m=l 

Notice also that is a subgroup of the twisted rotation symmetry group 50(4)' and 
that the lattice fields transform in reducible representations of this discrete group - for 
example the vector Aa decomposes into a four component vector and a scalar field 
(p under 50(4)'. Invariancc of the lattice theory with respect to these discrete rotations 
then guarantees that the theory will inherit full invariance under twisted rotations in the 
continuum limit. 

Proceeding in this manner it is possible to assign all the remaining fields to links on 
the A2 lattice. Since V'a(n) is a superpartner of Wa(n) it must also reside on the link 
connecting n ^ n + Ca- Conversely the field ul{n) resides on the oppositely oriented link 
from n — )■ n — e^. The ten fermions Xabi^) are then chosen to reside on new fermionic 
links n + + e„ — ^ n while the singlet fermionic field ?7(n) is assigned to the degenerate 
link consisting of a single site n. 
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The action of the theory takes the fohowing form 

^ n,a,b,c,d,e 

-^Tr CabcdeXdein + fl^ + p-h + /ic)^l^"''Xa6(n + Ac)} • (4-7) 
where the lattice field strength is given by 

TabH = --Vi+^Ubin) = --(l(a{n)Ub{n + fij -Ub{n)Ua{ii + Pb))- (4-8) 



9 9 

and the covariant difference operators appearing in this expression are given by 

p(+)/(n) = ZY,(n)/(n + /ij - /(n)Z^,(n), (4.9) 

pW/rf(n) = ^e(n)/,(n + - fd{^)UM + Pd), (4-10) 
pt(-)/^(n) = /e(n)ZYt(n)_^t(„_^j/^(„_^j, (4.11) 

pt(-)/„,(n) = /,fe(n)ZYt(n - /ij - U^{n + /i^ + - /2J/,fe(n - /ij. (4.12) 

Notice that these definitions reduce to the usual adjoint covariant derivative in the naive 
continuum limit corresponding to Ua = I + Aa{x) + . . . and furthermore guarantee that 
the resultant discrete expressions transform covariantly under lattice gauge transformation. 
Furthermore, this use of forward and backward difference operators guarantees that the 
solutions of the theory map one-to-one with the solutions of the continuum theory and 
hence fermion doubling problems are evaded j27]. Indeed, by introducing a lattice with 
half the lattice spacing one can map this Dirac-Kahler fermion action into the action for 
staggered fermions [28|. 

It is important to realize that the vectors n, n + etc appearing in this action do 
not correspond to the positions in spacetime of sites and links of the original A'l lattice 
itself - instead they span an abstract hypercubic lattice whose sites and links are given by 
integer valued lattice vectors. (Which are related to the r-charges defined in the orbifold 
formulation |jl^.) These 4- vectors /i„ are defined as 

Ail = (1,0,0,0) 
^2 = (0,1,0,0) 

/ig = (0,0,1,0) (4.13) 

ii4 = (0,0,0,1) 

/is = (-1,-1,-1,-1) 

The integer-valued lattice site n can be related to the physical location in spacetime using 
the Al basis vectors e^. 

4 4 

R = a'^dj.u ■ n)e„ = a'^n^e^ , (4.14) 

u=l v=l 
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where a is the lattice spacing. On using the fact that X^^^m ~ 0, we can show that a 
small lattice displacement of the form dn = /x^ corresponds to a spacetime translation by 

4 4 

dR = a^(/ii, • dn)e^ = a'^ip,^ ■ tim)^u = aSm ■ (4-15) 

i/=i i/=i 

The supersymmetry transformations on the lattice fields are almost identical to their 
continuum counterparts: 



QUa{n 
QV'a(n 

QXa6(n 
Qry(n 
Qd(n 



^a(n) 





i(p(+)W,)t(n) 




(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 



After the Q-variation, as performed in the continuum, and integrating out the auxiliary 
field d, the final lattice action is 



S 



1 



9' 



Tr \{vi+H{n)] ' (vi-'Hin)) + ^ (pt{-)z^,(n)) ^ - Xa6(n)Pf;V5] (n 



n,a,6,c,d,e 



-r?(n)P^( Va(n) - -eabcdeXdei^ + 'p,a + 'P'b + Mc)^l^ ^Xafe(n + A^c 



(4.22) 



To see that this action targets the continuum twisted theory one needs to expand Ua about 
the unit matrix [IToll^ 



Ua{n) = -lN + iAa{n) 
a 

Uiin) = -lN-iAa{n) 
a 



(4.23) 
(4.24) 



While the supersymmetric invariance of the Q-exact term is manifest in the lattice theory it 
is not immediately clear that the Q-closed term remains supersymmetric after discretiza- 
tion. Remarkably, this can be shown using a remarkable property of the discrete field 
strength which can be shown to satisfy an exact Bianchi identity just as for the continuum 

Hi- 

eabcde'Dl^-^J'Un + nc) = (4.25) 



5. Renormalization - General Analysis 

Power counting reveals that the continuum four-dimensional theory has an infinite number 
of superficially divergent Feynman diagrams occurring at all orders of perturbation theory. 

To leading order this is equivalent to the more conventional expression Ua{x) = ie""-^-'"). We will see 
that the linear representation offers important advantages over the exponential in our later calculations. 
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Of course in the continuum target theory all of these potential divergences cancel between 
diagrams to render the quantum theory finite. However, since the lattice theory does not 
possess all the supersymmetries of the continuum theory, it is not clear how many of these 
will continue to cancel in the lattice theory. 

Before we embark on a general perturbative analysis of this lattice theory it is instruc- 
tive to try to ascertain what kinds of counter terms are permitted by the lattice symmetries. 
In the case of lattice, these symmetries are 

a) Exact Q supersymmetry. 

b) Gauge invariance 

c) 1S5 point group symmetry and discrete translations. 

In fact, other than exact lattice supersymmetry, the U {N) lattice gauge theory also has a 
second fermionic symmetry, given by 

r/(n) r){n) + eIn, (5(all other fields) = (5.1) 

where e is an infinitesimal Grassmann parameter. Thus, we extend our list to include 

d) Fermionic shift symmetry 

In practice we are primarily interested in relevant or marginal operators; that is opera- 
tors whose mass dimension is less than or equal to four. We will see that the set of relevant 
counterterms in the lattice theory is rather short - the lattice symmetries, gauge invariance 
in particular, being extremely restrictive in comparison to the equivalent situation in the 
continuum. The argument starts by assigning canonical dimensions to the fields = 1) 
[^] = I and [Q] = ^ where ^ stands for any of the twisted fermion fields {)^,i^a,Xab)- 
Invariance under Q restricts the possible counterterms to be either of a Q-exact form, or 
of Q-closed form. There is only one Q-closed operator permitted by the lattice symmetries 
and it is already present in our bare lattice action. A possible renormalization of this 
fermion kinetic term is hence allowed. Beyond that the exact lattice supersymmetry forces 
us to look at the set of Q-exact counterterms. 

Any such counterterm must be of the form O = QTr ("^ f{U,U'^)). There are thus no 
terms permitted by symmetries with dimension less than two. In addition gauge invariance 
tells us that each term must correspond to the trace of a closed loop on the lattice. The 

smallest dimension gauge invariant operator is then just Q{Tiil^aUl). But this vanishes 

t • • • 7 

identically since both Ua and tpa are singlets under Q. No dimension ^ operators can 

be constructed with this structure and we are left with just dimension four counterterms. 
Notice, in particular that lattice symmetries permit no simple fermion bi-linear mass terms. 
However, gauge invariant fermion bi-linears with link field insertions are possible and their 
effect should be accounted for carefully. 

Possible dimension four operators are, schematically, 
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L2 = 5"'QTr ivT^iUa) 

Li = g-^QTr (r?)TV {UaU},) (5.2) 

The first operator can be simplified on account of the antisymmetry of Xab to simply 
QiXab^ab)-, which again is nothing but one of the continuum Q-exact terms present in 
the bare action. The second term in (|5.2| ) also corresponds to one of the Q-exact terms 
in the bare action. However the third term L3 is a new operator not present in the 
bare Lagrangian and the same is true for the final double-trace operator L4. Both of 
these operators transform non-trivially under the fermionic shift symmetry, but a linear 
combination of the two 

D = Li- (5.3) 

is invariant under the shift symmetry with the rank of the gauge group U {N) . 

By these arguments it appears that the only relevant counterterms correspond to 
renormalizations of operators already present in the bare action together with D. This 
is quite remarkable. The most general form for the renormalized lattice Lagrangian is 
hence 



n,a,b,c,d,e 
CK4r 

Y 



C= Yl { 2 Tr [ - iaiXabVi+^Kbin) - ia2vH'Di'--^Ka{n) + ^r]{n)d{n) 

' ,c,d,e 

-Tr eabcdeXdein + p-a + flh + p,jVl^~^Xab{n + /ij} + QPD , (5.4) 



where {ai,i = 1 . . . 4) and /3 are dimensionless numbers taking values (1,1,1,1) and 
respectively in the classical lattice theory. Thus it appears that at most four dimensionless 
ratios of these couplings might need to be tuned to approach = 4 Yang-Mills in the 
continuum limit. Furthermore, since these operators are dimension four we expect this 
tuning to be at worst logarithmic in the cut-off. 

In order to see the explicit form of the D operator close to the continuum limit, we 
expand the action around Z^m(n) = ^I. The result is 



L4 ~ — 

a 



5 5 
TV r?(n)(X] ^™(n)) - ^Tr 77(n)TV (^ V''"(n)) + 

m=l m=l 



(5.5) 



where ellipsis are dictated by super symmetry. The reader will immediately realize that 
(Sa=i ''l^a) is nothing but the ^5 (and twisted 50(4)') singlet contained in the reducible 
representation tpa- Indeed, it should be clear that it is the only field that could form a 
fermion mass term by pairing with rj. 

This is about as far as we can go by just using the lattice symmetries. We now turn 
to a full perturbative analysis to determine how the couplings (aj,/3) evolve with cut-off. 



6. Propagators and Vertices 

In this section we derive the propagators and vertices of the gauge-fixed A/" = 4, d = 4 
SYM theory on lattice. Then we write down the one loop diagrams relevant for the 
renormalization of the theory. 
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Upon rewriting the field strength and covariant derivatives in terms of the bosonic hnk 
fields Wa(n), the classical lattice action ( 4.22| ) takes the form 



S = Sb + Sp + Sc 



(6.1) 



where 

SB = \^Tr [(p(+)Z^,(n))^(p(+)W,(n)) + i (pt(-)z^^(„) 

n,a,b 

1 



-2 ^ TV (Ulin + naPlin) - Ulin + /i,)Z^,t(n)) (Ua{n)Ub{n + flj " l^binWn + fib) 



Sf 



and 



^ ^ ^ ^(^ac4d - Kdhc) Xab{n)(Uc{T^)lpd{T^ + Ac) " V'd(n)^c(n + (1^) 



5" ^ 2 



+^?(n)(Va(n)Z^i(n) -ZYi(n-/i„)V'a(n-/i„)) , 



(6.2) 



n,a,b,c,d,e 
X 



[xa6(n + /ic)^l(n) - ^](n + + M6)Xa6(n) 



(6.3) 



To proceed further we expand the Ua fields around unity as in eqn. [4.23| . Notice that this 
expansion point is but one of an infinite number of classical vacuum solutions - the full 
moduli space of the lattice theory corresponds to the set of all bosonic field variables Ua{n) 
such that 

= ^ Tr [(^,^(n + naMin) - Ul{n + fibMin)) {Ua{rv)Ub{rv + Ma) " Ub{n)Ua{rv + fib) 

n,a,b 

+ ^(l(ain)Ulin) -Klin- fiaWn - na))' 

These equations possess a large class of solutions corresponding to constant diagonal ma- 
trices modulo gauge transformations. We will use this additional freedom later when we 
compute the one loop contribution to the effective action of the theory. 



6.1 The Bosonic Propagators 

As usual it is easiest to compute the Feynman diagrams in momentum space. On the 
lattice a generic field ^'(x) has Fourier expansion 



$fxl = — - — \^ e^P''^. 



(6.4) 
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where x = a X]a=i ^a^a denotes the position on and the momenta he on the dual lattice 
given by p = Y^1=i "i^ada (fo^ ^ lattice with spacing a and length L). The dual basis 
vectors a = 1 ... 4 satisfy 

ea-db = ^ab (6.5) 

On an lattice both sets of lattice coordinates n^, rua take integer values in the range 
—L/2 + 1, . . . ,L/2. We will assume periodic boundary conditions in all directions in this 
paper. Eqn. 6.4 implies that fields are automatically invariant under translations by a 
lattice length in any direction and a field shifted by one of the basis vectors can be expressed 



as 



(6.6) 



where pa = ^rna- The only remaining is the question of how to deal with shifts in the lat- 
tice action associated with the additional 65 vector. However, the solution is simple: since 
Yla=i = we simply replace any 65 shift encountered in the action by the equivalent shift 
— Yla=i ^a- might have worried about an apparent lack of rotational invariance asso- 
ciated with the naive continuum limit of terms in the action which resemble Yla=i P<^ 
However, putting pa = p.e^ and taking the naive continuum limit this becomes 



(6.7) 



a=l 



^,1/ a=l 



which has the correct rotationally invariant form since the Greek indices refer to a Cartesian 
basis. 

Using these ideas the bosonic action when expanded around (4.23) and (4.24) gives 
the following second-order term in Fourier space 



« 2 ^ Tr (X(k) 6abfc{k)rAk) - f:{k)Mk) A(-k) 



k,a,b 



+B,{k) [/:(k)/fe(k)] B,{-k)) , 



where 



/a(k) = (e 



1). 



(6. 



(6.9) 



We need to gauge-fix the bosonic action before we derive the propagators. A natural 

.10) 



gauge-fixing choice would be an obvious generalization of Lorentz gauge-fixing |26| 

G(n) = ^(5(-)X(n)+a(-)X(n) 

a 

This gauge-fixing choice adds the following term to the bosonic action at quadratic order 

= ]^E^'H = ^ETr (5^)A(n))^ (6.11) 



n,a 



^For simplicity we will adopt the convention that momentum sums automatically include the 1 / {LaY 
normalization factor. 



-13- 



Figure 1: The bosonic propagator. 



wheiedi ^/(n) = /(n)-/(n-/i„). On using the relation X;n(^a^''/(n))5(n) = - En /H^a ''5'(n), 
the gauge-fixing term becomes 

Sgf = -- Tr Aa{n)di+^dt^A,{ii). 



(+). 



a 



n,a,b 



In momentum space it becomes 

5gf = - V Tr ^„(k)/:(k)/,(k)A(-k). 

rv ^ — ^ 



(6.12) 



5.13) 



Thus the gauge-fixed bosonic action to quadratic order is 



s'b^+Sgf « 2 Tr (Aa{k) J,fc/,(k)/*(k) - (l - -^)/*(k)A(k) A,{-k) 



k,a,6,c 



2a 



-fB„(k) [5,fe/,(k)/,*(k)] Bb{-k) 
The choice a = 1/2 makes the above expression diagonal 

Sg) « 2 J] TrX(k) [(5a5/c(k)/:(k)] A(-k) 

k,a,b,c 

= 2 ^ Tr [A(k),5„5(4 J^sin^ (|))A(-k) 



(6.14) 



k,a.6 



(6.15) 



Putting in the trace (using the convention Tr (T^T^) = ^6ab) the quadratic bosonic 
action can be written as 



(6.16) 



k,a,6 



where M^'^(k) = k'^dab^AB, with k^ = 4E^sin^ Thus only the AA propagator is 

non-zero and it is given by (See figure |l[) 



(^^(-k)X(k)) = KbSAB 



k2 



(6.17) 



6.2 The Fermionic Propagators 

The fermionic part of the action is of the form 

Sf = -\Y1 (xafc^fa^Vfe] + f?^l^"Va + ^eabcdeXab'^l^^hde) (6.18) 
abcde 
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Explicitly, we have 

Sf = -\ E Tr (xa6(n)Pf;V6](n) + r?(n)pt(-)^^(n) 



n,a,b,c,d,e 



1 



(6.19) 



when expanded up to second order in the fields using ( [4.23D and ( |4.24D , it becomes 



4'^ « 4 E Tr Xa6(k) [ - f:imc + /,*(k)^ Vc(-k) + r?(k)/,(k)V^,(-k) 

^ k,a,fe,c,(i,e 

+ ^ea6cdeXde(k)e^('="+'=''Vc(k)Xa6(-k) 



.20) 



Upon restricting the sum and rescaling the field 2xab ~^ Xab the fermionic action becomes 



4'^ ^ ^ E Tr (xa6(k) [ - /:(k)5,e + A* (k)5acl ^'^^"^^ + r?(k)/,(k)V^,(-k) 



k,a<6;c,d<e 



1 



+ ^ea6cdeXde(k)e*('''+'='')/c(k)Xafe(-k) 

We can then write this action in the form of a matrix product 



.21) 



5S?'.lE(*«*(-X))G)(_M°, 



Af(k) 
k) 



*(k) 
^-f-k) 



-L^$(k)Alci,(k) 
^ k 



(6.22) 



where $ = (*(k), ^'(-k)) and = (??, V'l, • • • , V'5,Xi2, • • • ,Xi5, • • • ,X45) and M(k) is given 
in block matrix form 





(77 Xde) (k) I -/:(k) 





/fc(k) 








-/J(k)(5ef, + /*(k)(5db eghcdcqghfcCk.) 



Xgh 



where qgh = e^^^s+^*^) , Notice that M has the properties M'^(k) = — M*(k) = -M(-k). 
Using the property that Yla /^a = can square the matrix to obtain 



M2(k) = - J^le^^'" - = -45;sin2 (^)li6 = -k 



a=l 



Thus 



a=l 



k2 



.23) 



.24) 



and the inverse of the full fermion matrix is 



1 / 



M^(k) 



= - — 

k2 V M(k) 



(6.25) 
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r?^(-k) 



Xfc(k) 



k2 



X^6(-k) 



Figure 2: The fermionic propagators. 



Then we can write the quadratic part of the fermionic action as 



53<S.(k)Aly(k)*,.(k) 



E <ff (k)A^.,(k)$f (k)Tr (T^r^) 



45 



k ij,A,B 



= E <ff(k)M,(k)$f(k)5^^ 

^ k ij,A,B 



(6.26) 



where we have expanded the fermions as $ = ^^T^ and used TV (T^T^) = ^6ab- Thus 
we write the propagators as 



or alternatively 



(6.27) 
(6.28) 



Notice that by switching the fields (with some relabeling) we have 



(*f (-k)*f (k)) = -(*f (k)^'f (-k)) = -^mU\,)5ba = -Z-Mi^{\^)5AB . (6.29) 



k2 



For a consistency check we replace k with — k and get 



k2 



(-k)*f (k)) = ^Mlji-k)SAB = -^Mijik)5AB . 



k2 



k2 



(6.30) 



We must also undo the earlier rescaling of the % field giving a factor of ^ in the ipx 
propagators and a factor of ^ in the XX propagators. It is also important to note that if 
we switch the direction of fermion flow in the propagators then we pick up an additional 
minus sign. 



-16- 



6.3 The Vertices 

Before we write down the expressions for vertices, let us further fix our conventions on the 
trace algebra. For the generators of U (N) one has 

T^T^ = ^{dABC + ifABc)T^. (6.31) 

where (Iabc and Jabc are the symmetric and antisymmetric structure constants, respec- 
tively. This product formula is consistent with our previous trace convention TV (T^T^) = 
\6ab and in addition yields the results 

TV {T^T^T^) = Tr (^^{cIabd + ifABD)T^T^^ (6.32) 
= \{dABD + ifABD)Tr [T^T^] 

= -^{dABD + ifABD)-^SDC 

= -^{dABC + ifABc) = ^AabC ■ 

Since /abc is antisymmetric and dAsc is symmetric it follows that 

Aacb = Aasc ■ (6.33) 

To extract expressions for the vertices we now return to the original gauge-fixed action for 
the theory given by 

^ = i E ^ [K^Z^.(n)) ' {vi-^Hin)) + \ (Pi(-)Z^„(n)) " 



3 K W 

n,a,o,c,a,e 



+ 



2A„(n)9(+)9j-)A(n) - (xa6(n)P;;V(,] (n) + r?(n)pt(-)^„(n) 



+ ]^(^ahcdeXdeiP- + lj'a+lj'b + V'cWl^ \ab{T^ + Tl- S) ■ (6-34) 

The last three terms of the action give rise to vertices between varying number of ^'s and 
the fermions r/, i/'a; and Xab- There are three vertices that arise at linear order in A: 

• The %l;Ar] vertex 

n,o 

= - ^ TV (r?(n)^„(n)Wi(n) - r?(n)Wt(n - /iJVa(n - /ij) 

= _ ^ TV e*(''+'i+P)-"(,7(k)V'„(q)(-i)A(p)-??(k)(-i):4„(p)e'^'"Va(q)e*''") 

n,k,q,p,a 

5-k,q+p^^(k):4f (p)Vf (q) 5ab[>^ABC - \ABce-'^-+^^^] ■ (6.35) 



k,q,p 
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i^np) p 



-T(-^ac^6d + SadSbc)[><ABCe''''= - XABCe'''''] 



xtbip) P 



-k Xd^(k) 



o ^abcde \ ^ 



Figure 3: The vertices connecting fermions and complexified gauge fields. 

Thus the Feynman diagram contribution for this vertex is (add a minus since it comes 
from the first order term of e"*^) 



(6.36) 



• The tpAx vertex 

Vi'Ax = -J2^ Xa6(n)P[+V6](n) 



n,a,b 



= 5^ IV ( - x„6(n)p(+)V6(n) + Xa6(n)Pi+Va(n)) 

n,a,b 

XI i-^acSbd + Sad^bc)"^ Xabi^) (^c(n)V'd(n + Mc) " '^di^)Uc{n + Md)) 



n,a,b,c,d 



= ^-Ka+pi-^acSbd + SadSbc)Xabi^)A^i^)^diP)l[^ABCe'P^ - XABce"''] • (6.37) 

k, q,p, a, 6, c,(i 



The vertex is given by 



^X^V- = --^{-^achd + 5ad5bc)[XABCe'^'' - AABCe'*'^] . 



(6.38) 
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The x^X vertex 



V 



n,a,6,c,d,e 



= Yl ^abcdeXdei^ + fia + 'fib + K) (xa6(n + Mc)^c (») - U^ci^ + fia + fibJXab^) 



n,a,6,c,d,e 



E '^-k,q+pWexS.(k)-4f (q)x^b(p) ( e^(^"+'=^+'==)^[A^Bce^^^ - AABce^^^^+'^^^l 



\s.,p,ci,a.,b.,c,d,e 



The vertex is given by (taking into account both possible contractions with external 
propagators) 

yxAx = -y-bcde(^^'^''''^''''^'''^^^^BCe'P'-X^^ 

(6.39) 



7. One Loop Diagrams for the Renormalized Fermion Propagators 

Using these propagators and vertices it is straightforward to see that the renormalized 
fermion propagators receive contributions from the following four amputated diagrams. 

• The amputated ijip diagram. We have an AA propagator, a V'X propagator, an r]A^p 
vertex, and a x^''P vertex. Using the expressions above we have 



^.^P) = 2;2;i;5-p,k+q[i[(e''=^-l)<^ac-(e^'==-l)5„6] 

k,q BC abc ^ 

• [^(-56a<^cd + SMSca)[XBCDe-'^'' - XBCDe''^''] 



-[Xabc - XABce'^"''^'''^^] 



(7.1) 



• The first amputated ipx diagram. We have an AA propagator, a XX propagator, a 
ipAx vertex, and a X-^X vertex. 



4x(P) = E E E[iiWe^^^^"'^^^(e^^--l 



k,q bcdefm BC 

1 



' [ ~ 4^~^b^^'^ + SbaScd)[XACBe'^'^ - Xacbs 
■eefd9h(e"'('^+''+''^[XBCDe-'P'' - ABCDe'^^^+«'')] 



-iqa 



(7.2) 
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(i) "0 — 77 diagram 



(ii) First ?A — x diagram 



A^(-q) 



(iii) Second ijj — x diagram 
-Af(q) 




^^iv) r;^(-k) (k) x^e(-P) 



^!r(-q) 



(iv) X ~ X diagram 

^?(q) 




xi,{p) Xd^e(-k) ^f(k) x,^,(-p) 



Figure 4: The one loop diagrams of fermions and complexified gauge fields. 



The second amputated V'X diagram has an AA propagator, an rjijj propagator, a il^Arj 
vertex, and a ipAx vertex. This yields 





1 








■ Sab 









k,q be BC ^ 

■ [ - ^{-Sdb^ec + SdcSeb)[>^DCBe'^'' - XoCBe''^' 



[Aacb — Xacbc 



-i{Pa~qa)] 



(7.3) 



The amputated XX diagram. It has a AA propagator, a xV' propagator, a X-^X 
vertex, and a ipAx- 

^xx(p) = E E E -^k+q-P.o - - (e-*''^ - l)5fe 

k,q cde/ BC ^ 





1 







-i{ga+gb)] 



(7.4) 



In appendix A we show that the contributions of these diagrams all vanish in the limit 
p — )• indicating that mass counterterms are absent in the lattice theory at one loop. In 
our general argument of section |5| we argued that the only dangerous mass term involved 
a coupling of r] and -i/^a- We now see that this term does not arise at one loop. In the next 
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section we will show that this feature persists to all orders and thus our general conclusion 
will be that no mass counterterms are needed at any finite order of perturbation theory. 



8. Effective Action 

In this section we will compute the partition function of the lattice theory in one loop order 
around an arbitrary classical vacuum state in which the fermions vanish and the bosonic 
fields correspond to constant commuting matrices. To start we expand the fields around 
such a constant commuting background, 

l(a{n)=Ua + iAa{n), = - iAa{n) (8.1) 

Choosing the gauge a = 1/2, the quadratic part of the bosonic action then takes the form 

Sb = -2J2 Tr A(n)pt(-)p(+)^,(n) . (8.2) 



n,a,f) 



iere the covariant derivatives depend on the constant commuting classical background 
UaMa = 0. After integration over the fluctuations in the bosonic fields one finds the 
aosonic contribution to the one loop partition function is given by 



det-5(pt(-)p(+)) (8.3) 
The gauge fixing functional ( 6. 10] ) leads to the quadratic ghost action 

5G = 5]TrcPt(-)p{+)c. (8.4) 

n,a 

The quadratic fermionic part of the action is given by the corresponding terms in ( 4.22| ) , 
except that now the covariant derivatives depend only on the background fields. 

Since the background is constant, we can pass to momentum space in which the action 
separates into terms for each mode k. The 16 x 16 fermion matrix M(k) for the mode k 
then can be shown (using MAPLE to compute the determinant) to satisfy 

det M(k) = det(pt(-) (k)Pi+) (k))*^ . (8.5) 

Going back to position space, and taking into account the fact that there is a double 
counting of modes in the matrix form (|6.22|), we obtain 



Vi{M) = det^(pt(-)p(+)) . (8.6) 

The ghosts add another factor of det{T>]j ^V^^), which is just what is needed to cancel 
the bosonic contribution given earlier. 

In conclusion, we have shown that the one loop effective action of the lattice theory 
obtained by expanding about an arbitrary point in the classical moduli space is identically 
zero. Thus, as for the continuum, the moduli space is not lifted in this analysis and 
hence there can be no boson or fermion masses at one loop. Furthermore, we expect 
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that we can extend this analysis to all loops since the partition function of the lattice 
theory is a topological invariant and hence can be computed exactly in the semi-classical 
approximation (see Appendix C). Indeed, Matsuura uses similar arguments to show that 
the vacuum energy of supersymmetric lattice theories with four and eight supercharges 
remains zero to all orders in the coupling [30|. The calculation presented here extends 
this to the case of sixteen supercharges^. Thus we conclude that boson and scalar masses 
remain zero to all orders in the coupling constant. This implies that the fermions also 
remain massless which is consistent with our explicit one loop calculation. 

At this point we have derived expressions for the amputated one loop diagrams that 
contribute to the renormalization of the three twisted fermion propagators. This is sufficient 
to calculate ai, 02 and 04 that appear in the general action 



E {2Tr[ - iaiXabV^^^Ubin) - ia2r/(n)pt(-)^/^(n) + ^7?(n)(i(n) 

^ n.,a,b,c,d,e 

- yTr eabcdeXde{n + p.a + 'P'b + Mc)^J^"-'Xa6(n + Tic)] ■ (8-7) 

However the coefficient 03 requires further work. One simple way to extract it is via a 
computation of the renormalized auxiliary boson propagator which we turn to in the next 
section. 



9. One Loop Diagrams for the Auxiliary Field Propagator 

We have shown that the off-shell form of the bosonic action is given by 

5b = E {:Fl{n)Tab{n) - ^d(n)pt{-)zY„(„) + ^d\n)) , (9.1) 

n,a,b 

where Jab(n) = -pi+^Ubin). 

In our previous computation of the fermion diagrams we integrated out the field d to 
give an on shell action defined just in terms of the complex gauge link fields Ua and Ua- 
In this section we will not do this but instead focus on a computation of the renormalized 
propagator for the d field. The Feynman rules for the fermions will be identical to our 
previous scheme but the boson propagators will change and so we need to recompute those 
propagators in this off shell scheme. We proceed in the standard fashion by expanding the 
link field Z^a(n) 

Ua{n) = l + igAa{n), Ul{n) = 1 - igAaH- (9-2) 
and using the same lattice gauge-fixing term as before 

Sgf[A] = --^Tr (9(-)^,(n))2, (9.3) 

n,a 

''Notice that in this calculation we have not included any mass terms that would guarantee the stability 
of the initial classical vacuum state we have chosen to expand around. We have also ignored a potential sign 
problem associated with the replacement of a Pfaffian with a square root of a determinant. Nevertheless 
we expect the result to be robust; the existence of an exact supersymmetry should ensure that the object 
we are computing is a lattice regularized Witten index and hence independent of both coupling constant 
and background field. 
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we find the momentum space form 

Sgf[A] = - V Tr A(k)/:(k)/,(k)A(-k) . 

rv ^— • 



(9.4) 



It is convenient in tliis calculation to work with the real and imaginary parts of the complex 
gauge field explicitly. Thus 

Aa = Aa + iBa (9.5) 

The gauge-fixed bosonic action on the lattice to quadratic order in fields, with the choice 
a = |, is then 

Sf = ^ Tr 2^,(k) [<5,fe/e(k)/*(k)] A6(-k) + 2Ba{^) [5„,/,(k)/:(k) - /:(k)/6(k)] i?,(-k) 

k,a,f) 

-2id(k)/,(k)i?„(-k) + ^d(k)d(-k) (9.6) 
We see that the d — Ba system decouples from Aa to this order. Its action is given by 
S^^\d,Ba] ~ J]Tr2i3,(k)[Wc*(k)/c(k)-/:(k)/fc(k)ji?fc(-k) 



k,a,b 



-2id(k)/,(k)i?„(-k) + -(i(k)d(-k) 



or in matrix form 



[d Ba) (k) 



-i/*(k) M,fe(k) 



d 

Bb 



(-k) 



(9.7) 



(9.8) 



where Mafc(k) = 2[(5ab Ylc /c(k)/*(k) — /*(k)/f,(k)]. Using standard identities for the inverse 
of a partitioned matrix we find 



I -iA(k) 
-i/*(k) M„fe(k) 



1 



i/6(k) 



E,/c(k)/*(k) I i/*(k) 



2-^5 



(9.9) 



Thus the lattice propagators are 

(d^(k)d^(-k)) = 
(d^(k)Sf (-k)) = i^AB- 



We have Ec/c(k)/c(k) = 4^^sin2 and as before we define k^ = 4^^sin2 



-ika _ J") 



k2 



2k2 



From eqn |9.q the propagator for the A field is also 



{A^{k)Aj^i-k))=6ab5AB 



2k2 



(9.10) 
(9.11) 

(9.12) 
(9.13) 



-23- 




Figure 5: The generic diagram contributing to the renormalized d propagator 



Notice that the field d is non-propagating at tree level. Using these propagators and those 
derived earlier for the fermions and ghosts we can now write down the generic Feynman 
diagram contributing to a renormalization of the auxiliary boson propagator. It is shown 
in figure |5| and represents the set of amputated diagrams possessing two external B field 
legs. These combine with the external < dB > propagators derived above to yield the 
renormalized propagator for the auxiliary field d. Notice that the vanishing of the tree 
level < dd > propagators ensures that no amputated diagrams with 2 d field external legs 
contribute. The set of all such lattice Feynman diagrams is shown below and corresponds 
to a subset of the B field vacuum polarization diagrams. It is important to notice that 
almost all these diagrams appear in the continuum off shell twisted theory; the exceptions 
are just the diagrams containing a BBd vertex which corresponds to the lattice vertex 

VdBB = (d^(-k-q)i?f(k)5f(q)) = l5„,(A^BC + W)(l-e-*('='"+''»)) • (9.14) 

Clearly this vertex vanishes as the lattice spacing is sent to zero and hence this diagram 
does not contribute to the divergent piece in the < dd > propagator at this order of 
perturbation theory. 

Hence we are left with a set of diagrams which correspond to those of the equivalent 
continuum theory at one loop order. This fact can be exploited later to allow us to argue 
that the leading logarithmic divergences of the lattice theory are shared with the continuum 
theory. Anticipating this we will not write down explicit expressions for these amputated 
lattice diagrams in this section. 



10. Divergence Structure of the One Loop Diagrams 

At this point we have derived expressions for the amputated one loop diagrams that deter- 
mine the renormalization of three fermion propagators and also the set of Feyman graphs 
needed to renormalize the auxiliary bosonic field propagator. In principle, this input will 



allow us to determine all four coefficients a, appearing in the renormalized action eqn. 5.4 . 
Of course the question of how much fine tuning is required to regain full supersymmetry is 
determined by the parts of these expressions which diverge as the lattice spacing is sent to 
zero. We must therefore evaluate the expressions for the one loop integrals as the lattice 
spacing tends to zero. 

First, let us discuss the diagrams contributing to the fermion propagators. We have 
shown in Appendix B that the three fermion amplitudes all vanish for vanishing external 
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Sc^(-q) 



Sa^(p) 



Bf(k) 




i3f(-q) 




i3f(q) 



sf(-p) 



A^(-k) 



c^(-q) 



Sa^(p) 



c^(k) 



c^(q) 



S.f(-P) 



c^(-k) 



B^ip) 



*^(k) 



^f(-q) 




Sa^(p) 



*^(q) 



BFi-p) 



l'^(-k) 



^f(q) 



sf(-p) 



i3f(-q) 




sf(q) 



5a^(p) 



5f(-p) 



B^ip) 



^f(k) 




i??(q) 



i3i^(-p) 



^f(-k) 



i3^(p) 



Sf(k) 




S.^(-P) 



rf^(-q).- 



Bf(q) 



B^ip) 



i?f(k) V 



S.r(-P) 
rf^(-k) 



Figure 6: Set of all lattice amputetated Feynman diagrams contributing the renormalized d 
propagator 



momentum which is consistent with our effective action computation showing that no 



fermionic mass terms can be generated perturbatively. Reisz's power counting theorem |31]- 
|34] shows us that we cannot simply take the naive continuum limit of the expressions for 
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the amputated one loop diagrams as they have a naive degree of divergence of 1. However 
we can use a trick due to and detailed in to extract the leading divergences. 
We split the integral I(p) into two pieces as follows: 



lim/(p) 

a— 5-0 



lim 

a->-0 



/(p) - /(O) - ^ 



Pb 



di_ 

dpb 





+ lim 


p=0 


a-s>0 



^ OPb 



(10.1) 



The first term in square brackets can now be evaluated in the naive continuum limit and 
contains no divergence. The second term contains the divergence but contains no external 
momenta in the integrand which simplifies its evaluation on the lattice. In addition we 
know that /(O) vanishes for each of our diagrams so the calculation becomes simpler still. 

We will find that the resulting expressions have logarithmic divergences of the form 
In ij,a where is a small mass parameter used to regulate the behavior of the integrand 
close to the origin of momentum space and a the lattice spacing''. 

One obvious way to proceed is simply to numerically evaluate the integral for a variety 
of regular masses fj, and extract the logarithmic divergence and any constant contributions 
using a fitting procedure. However, if we are only interested in the leading log divergences 
there is a simpler approach detailed in the next section in which a naive continuum limit 
can be taken and the expressions evaluated using, for example, dimensional regularization. 

In the next section we give an example of this procedure for the amputated rjip diagram 
and show how to extract similar results for the remaining fermion self energy diagrams. 
We will also see that the same procedure allows us argue that the leading log divergent 
contribution to 03 is also equal to its value in the continuum theory. 

10.1 The Amputated Fermion Diagrams 

We start with our simplified expression for I^^^ (p) given in Appendix A 



^vi^d (P) 



(2^ 



E 

BC 



_e*(P-'?)'^) 

(p - q) 



E [dABcdscDie-'P^ - e"^'' - 1 + e'P^+"''') 



+fABcfBCDie~ 



■^Pa _^ g^gd + 1 + e 



ipa+lQd 



)] 



(10.2) 



As a first step we need to calculate the derivative of the diagram (re-inserting the lattice 
spacing a and the infra-red cutoff fi) 



dlri^bd (P) 



dpb 



p=0 



^ (27r)4 (q2 + ^2a2' 



o. d^q 



-a 



(27r)4 (q2 + ^2^2)2' 



^ (1 _ e— 9d)/^Bc/i?CD(l + e*"^'^) 



-ia5dbe-"''''')fABcfBCD{l + e^'^'^'^) 



- j4 

a a q 



-a 



[l-e 



(27r)4 8(q2 + ^2^2)2 
X '^{dABcdsCD + fABcfBCD)Sabi-ia){l - € 



laqa ^ 



(10.3) 



^We will only consider the case of infinite lattice size which reduces all lattice sums in momentum space 
to integrals. 
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A further simphfication now occurs; if we are only interested in the leading log fia coef- 
ficient we can evaluate this integral in a small q region around zero. This is because the 
contribution of the integrand to the log fia coefficient comes only from small q. Further- 
more in the region q the propagators and vertices inside the integral will approach 
their continuum counterparts and hence the logarithmic divergence can be extracted by 
replacing the lattice integrals by their naive continuum limit. Note that this only works 
for the coefficient of the log - we must evaluate the integral numerically (and then fit) in 
order to extract the constant terms. This (longer) calculation is in progress and will be 
published in a future paper. Thus we find: 



lim 



oo j4 



d^q -4:iqbqd 
p=0 ■ ^-oo (2vr)4 {q'+n'r 



fABcfBCD 



f°° d^q 2i 

+ / T^TUTT', 2^^dbfABcfBCD ■ (10.4) 



Note that we cannot just set the first term in this expression to zero as and Cb are not 
orthogonal to each other, instead we have 



/ 



Then Cb ■ Cd = S^d — ^- Wc use dimensional regularization and the fact that Ylb P& = to 
evaluate the resulting integrals getting 



Pb 



dlv^^AP) 



6 



p=0 



~ --^PdfABcfBCD^ogna . (10.6) 

Note that we have inserted the cutoff - inside the logarithm to ensure that it is dimension- 
less. 

Since all the Feynman graphs we need to evaluate are logarithmically divergent and in 
one-to-one correspondence with continuum diagrams, the resulting logarithmic divergences 
can all be extracted by following a similar procedure i.e. taking the naive continuum limit 
of the relevant /(p). 



^^O^'iLM - j-^^Y. ^^TJmP^WTW) ^^^^'^"^^ " ^^"^^""'^ f ABC f BCD 

3i 



-^^fABcfBCD {SagPh - ^ahPg) log lia (10.7) 



(2) ^ ^ f cg^q -ijjp - q)c 

7r)4 4-2(g2 + ^2)((p_q 

f ABC f BCD (SdaPe - ^eaPd) log Ha (10. 



lTo^V^"X.e(P) - I E 2(g2 + ^2)'^p _'q)2 + ^2) iS<ia6ec " Sd.6,a) IaBC IbCD 



327r^ 
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► C ^ 5 ^ ^ C ^ L ^ 

Ipd V i'c V i'd Xgh V'c 



Figure 7: Full rjip Propagators 



This obviously leads us to define I^„y. (p) = \, (p) + /i (p) and therefore 

1 

^i'aXdeiP) ~ -^fABcfBCoiSdaPe - ^eaPd) 



(10.9) 



lim /i 




2((?2+/i2)((p_q)2 + ^2) 



GabdghfABcfBCD 



-{h^g) 



f ABC f BCD 2. (^abdghPd fog 



(10.10) 



Note that these calculations of the log terms for the other diagrams have also been verified 
by numerical evaluation and fitting of the resulting lattice integrals. 

10.2 The Auxiliary Field Diagram 

Since the amputated divergent diagrams for the lattice d propagator are log divergent we 
can extract the sum of these logarithmic divergences using the same tricks we used for 
the fermions; evaluating the diagram in the naive continuum limit. The sum of all these 
diagrams, contracted with external dB propagators, will then yield a log divergent term of 
the form 



where c is a constant to be determined by explicitly evaluating the diagrams. However, 
we will argue in the next section that there that it is not necessary to evaluate these 
diagrams, even in the continuum, to determine 03 - the requirement that the continuum 
theory preserve full supersymmetry will automatically determine in terms of the other 
ai corresponding to the fermion propagator renormalization. 

10.3 Prom Amputated Diagrams to Renormalized Propagators 

The leading logarithmic divergences appearing in the renormalized propagators are ob- 
tained by combining the (divergent parts of) the individual amputated diagrams we have 
just computed. In principle several of the amputated fermion diagrams can appear as in- 
ternal bubbles when correcting a given fermion propagator. As an example consider the 
iprj diagram shown in Figure 0. Naively we see that three of our amputated diagrams con- 
tribute to the renormalization of this propagator. However we find that (at least in the case 
of the log divergences) the Lorentz structure of the propagators and integrals means that 
only the "qip amputated diagram contributes to the renormalization of the r]ip propagator. 



Cdd = cJacbIdcb In (/ia) 
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We demonstrate this through explicit calculation. Denoting the full diagrams by C 
and noting that as we are dealing with only the divergent part we can approximate the 
lattice propagators by their continuum analogues we find 



Wd*? - ^ Z^-''?V'clPj— 2- + 2^ -2 ^XshV-clPj— r 

c,g,h 



-^^fABcfBCD 2^ Pc^ log ixa 

~ -^fABcfBCD^^ log na . (10.12) 
The second term disappears as 

'}2{PgSdh - PhSdg){ScgPh - SchPg)Pc = ^{PcPd - ^cd - P^^dc + PcPd)Pc 

c,g,h c 

= . (10.13) 

We can similarly show that only /^^ contributes to C^^^ and I^^ to C^^. Note, however, 
that this analysis strictly only applies to the logarithmically divergent piece in C. 

^ iff iPgSah - iPhKg ^ X X JPd^ce - iPe^cd , 
'^^aXde = ^JABcfBCD ^2 ^^QcPh - dhcPg) ^2 

g,h,c 

1 . f iPd^ae - iPeSad , /1n1^^ 
= -^fABcfBCD ^2 logyua ■ (10.14) 

In calculating C^-y- we must take into account that the internal propagator in /^^ can be a 
V'X or xV'- This contributes another factor of 2 to C^^. 

^XabXde = --^JABCjBCD^OgUa ^abcfg-^^fgihjPi^hjkde-^^ 

c,f,g,i,h,j,k 

= --^fABcfBCD log/xa^ '^P^Y^ {^ai^bckde + ^U^cakde + ^ci^abkde) 

c,i,k 

= f ABC f BCD log fia'^^eabkde ■ (10.15) 

The coefficients are now determined by the coefficient of the propagator in the renor- 
malized propagator amplitudes C. Explicitly we find 

ai = l + 6jln/xa z = 1,2,4 (10.16) 

where bi = b = Note that we have used fABcfscD = ^^ad- This is required as 

the colour structure of any counterterms must match the tree propagators. However this 
is strictly only true for SU{N) as fABcfBCD = N{6ad — (^ao^do) for U{N). This does 
not matter in the continuum as the U{1) trace piece simply decouples from the rest of the 
system and can be ignored. When doing lattice simulations we might imagine achieving 
a similar result by giving the U (1) mode a large mass of the order of the cut-off which 
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will serve to decouple it from the SU{N) modes at finite lattice spacing. The breaking of 
supersymmetry in this sector may then be removed by sending this U (1) mass to zero after 
taking the continuum limit. 

While naively one might have expected the coefficients bi to be all different our re- 
sults indicate that in fact the log divergent parts of bi and hence ai are actually all equal. 
This fact can be understood quite simply; to untwist the continuum theory into a the- 
ory with four Majorana spinors requires that the continuum twisted fermions exhibit a 
common wavefunction renormalization. This just follows from the fact that the individual 
components of the spinors mix the different twisted fermions together. To achieve this 
requires that the corresponding renormalization constants of the kinetic terms Oj should 
all be equal - just as we find. Furthermore, since the leading log behavior of the lattice 
theory is the same as the continuum we should expect that the log divergent part of the 
lattice couplings behave in the same way. Thus a single wavefunction renormalization of 
the twisted lattice fermions is all that is needed to render the renormalized theory finite. 
The common anomalous dimension of the fermions in this twisted scheme is then given by 

In the case of the < dd > propagator the leading log divergent contribution can be 
computed from the naive continuum limit of the corresponding continuum expression for the 
sum of the BB bubble diagrams given in diagram ^. Combined with the fact that the tree 
level < dB > propagators required on the outside of these BB amputated diagrams are the 
same as the continuum to 0{a) we find that the log divergence in the mass renormalisation 
of the d field must be the same on the lattice as in the continuum. Using this fact we can 
argue that the log divergent part of 03 must actually be equal to that of the fermions eg 
ai. This follows from the fact that the bosonic action for general ai can be rewritten as 

ai(TabJ'ab)+^^[l[Da,Vaf) (10.17) 

Only for = 02 = ai can this renormalized bosonic action be untwisted to yield the 
conventional gauge field plus scalar action in the continuum limit. But since the continuum 
twisted theory possesses full supersymmetry this must be true. And our general arguments 
then tell us the log divergence of as on the lattice must satisfy the same property. 

To summarise; we find that the log divergent parts of the coefficients ai,i = 1 ... 4 
must all be equal to one loop order in the lattice theory. This implies that a common 
wavefunction renormalization of both twisted fermions and bosons is sufficient to render 
the renormalized theory finite at one loop with all fields acquiring an anomalous dimension 
(in this scheme) given by 7 = Physically, the equality of the couplings ai,i = 1 ... 4 
means that no logarithmic fine tuning is required at weak coupling for the lattice theory 
to exhibit full supersymmetry as the lattice spacing is sent to zero. 

11. Conclusion and Discussion 

We have examined a recently developed lattice construction for A/" = 4, d = 4 super 
Yang-Mills theory using perturbation theory. We argue that the exact symmetries of the 
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classical lattice theory; namely gauge invariance, a single exact supersymmetry Q and the 
(large) point group symmetry of the lattice strongly constrain the possible counter terms 
induced by quantum corrections. Indeed with one exception the only relevant counterterms 
correspond to renormalizations of existing terms in the action. We furthermore show by a 
computation of the effective action that the one new operator which cannot be excluded 
in the general analysis actually makes no appearance to all orders in perturbation theory^ 

The renormalized action can then be written in terms of 4 coupling constants ai which 
take the value unity in the classical lattice action. We evaluate the renormalization of 
these couplings at one loop using lattice perturbation theory. Three of the couplings can 
be computed by examining the renormalization of the three twisted fermion propagators. 
The final coupling is most easily read off from a one loop contribution to the propagator 
for a bosonic auxiliary field. The relevant propagators and vertices arc derived and the 
amputated one loop diagrams constructed. All these diagrams possess identical logarithmic 
divergences of the form In /xa where a is the lattice spacing and /j, a mass scale introduced to 
regulate the small momentum behavior of the integrands. This divergence can be absorbed 
by a common wavefunction renormalization Z of the twisted fermions and bosons. 

The simplest way to understand this rather surprising result is to realize that the 
coefficient of the logarithmic divergence of some one loop diagram in the lattice theory can 
be extracted by taking a naiVc continuum limit of the diagram, since the log divergence 
comes from the small loop momentum region of the integral. Provided that the lattice 
diagrams correspond one-to-one with equivalent continuum diagrams,^ and that all lattice 
propagators and vertices reduce to their continuum counterparts for small momenta, this 
means that the log divergences of the lattice theory are equal to the same divergences in 
the continuum theory^'^. Furthermore, since the twisted continuum theory is equivalent to 
the usual = 4 theory in flat space it must possess the full Q = 16 supersymmetry. This 
fact ensures that all divergences present in the twisted continuum fermion self energies 
must be equal - which is indeed what we find. And this structure is necessarily inherited 
by the log divergent parts of the lattice theory at one loop. This is what leads to our main 
result; that only a one time tuning of the finite parts of the wavefunction renormalization 
needs to be performed at one loop in order to restore the full supersymmetry. 

This similarity between the divergence structure of the lattice theory and the con- 
tinuum theory is strongly suggestive that the beta function of the lattice theory will also 

*The calculation we have done is somewhat formal as it ignores possible instabilities associated with 
the flat directions and specifically the U{1) trace mode of the scalars. It is possible that regulating these 
directions by eg introducing a mass term for the (7(1) mode might modify our conclusions since supersym- 
metry is broken by such terms. Notice though that at large N the dangerous double trace mass operator is 
suppressed and hence the result should certainly hold in that limit 

®Notc that this is true in the present, twisted construction and that it is not true in other, more naive 
constructions. There, one would have 0{a) multigluon vertices that correct the fermion self energy through 
a tadpole. This lattice diagram does not occur in the perturbation theory we have described above. 

^"Clearly both continuum and lattice perturbation theories must employ the same scheme for this to be 
true. In this light we would note that it is possible by addition of sufficient auxiliary fields to arrange for 
Z = 1 for all fields in continuum J\f = 4 SYM. However such a superfield approach is not possible for our 
twisted construction. This does not spoil the finiteness of the theory which depends only on the vanishing 
of the beta function which is true in both twisted and supergraph schemes. 
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vanish at weak coupling. We leave a proof of this (and a computation of the finite parts 
of the diagrams) to a later paper and here merely give heuristic arguments as to why this 
result may indeed hold. First, note that the calculation of the beta function requires the 
evaluation of one loop vertex diagrams in the lattice theory. Preliminary calculations sug- 
gest that the set of relevant lattice vertex diagrams correspond one to one to continuum 
vertex diagrams and remain only logarithmically divergent. They may thus be evaluated 
in the continuum theory. The coefficient of this log divergence is then combined with the 
wavefunction renormalizations determined above to yield the one loop beta function in the 
usual manner. However we already know the result of this computation for the continuum 
theory; the beta function vanishes. We hence expect a similar result to hold at one loop 
in the lattice theory. Thus for weak coupling wc expect the lattice theory to possess a 
line of fixed points parametrized by the bare coupling constant just as for the continuum 
theory. However, our calculations do not reveal whether this feature survives in the lattice 
theory to strong coupling. At two or more loops the divergences of the lattice Feynman 
diagrams will not be equal to the those of the continuum theory and hence we cannot use 
the latter to infer the divergence structure of the lattice theory. To understand how to take 
the continuum limit in this regime will then require a mixture of two loop and numerical 
calculations. 
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A. Simplification of the One Loop Diagrams 

This section shows how the diagrams in the paper are simplified. For the easy evaluation 
of the diagram we use the following identities. 



^ABC^BCD = ^2 
B,C B,C 



d ABC d BCD — f ABC f BCD ■ 



(A.1) 



^ACB ■ 



(A.2) 



These relations imply 



{AX ABC — BXabc){CXbcd — DXbcd) 



dABcdBcniAC - AD - BC + BD) 
-fABcfBCD{AC + AD + BC + BD) . (A.3) 
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Starting with /,,^(p): 



(27r)^ 



EE 



16(p - q) q2 



((1 - e'^-''^'')5,a{SbaScd - SbdSca) - (1 - e'^^-'^^^){dbaScd - SMSca)Sba) 

7 (27r)4^^ 



a,b BC 



8(p - q) q2 



That is 



a BC 



[dABcdBCDie-'P'' - e'"^ - 1 + e'f<'+'«<^) 



_8(p-q) q2_ 

+fABcfBCD{e-'P'' + e^^-^ + 1 + e'P'^+"i'^)] (1 - e^(^'-^)")<5„d 
- [dABcdBCDie-'P-^ - e"''^ - 1 + e'P^+'i'^) 

+fABcfBCD{e-'P'' + e"i^ + 1 + e^^'»+^«<^)] (1 - e^(f-«)<^) 



y (2vr)^ 1^ 



^2 
8(p - q) q2 



i(p-q)d^ 



[dABcdBcnie-'P'^ - e''^'^ - 1 + e^(f+«)'') 



+fABcfBCD{e-'P'' + e^i" + 1 + 

- J2 [dABcdBCD{e-'P'' - e"i^ - 1 + e^f<'+^«'^) 

a 

+fABcfBCD{e-'P'' + e"i^ + 1 + e'P'^+"''')] 
1 



7 (2vr)^ 1^ 



_ p«(p-9)<i 



(,p - qj 



(A.4) 



Now / 



(1) . 



2 ^hcmeftghdefi ^ ' 



b,c,d,e,m,f B,C 64q2(p - q) 
{e'^"'^"^-l){SbdSca-SbaScd) 

X (A^Bce*^'* - AABce-*^") (e-'P('^+^+n (xBCDe'''^^+f^ - XBCDe'^"-"^'') 
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That is 



^ ^ rf,e,m,/B,C 16q2(p - q) 

^^""^ d,mB,C7 8q2(p-q) 
X (Sag^mh + ^ah^md^dg + ^ad^mg^dh — ^ah^mg — ^ag^md^dh — ^ad^mh^dg) 

X {dABcdBCDie''^'^~^''<-o+>'^^ - 1 - e'i<-o+>^-) + e-^(P<^+«")) 
-fABcfBCD{e'^''^'i^<'+'-)^ + 1 + e''?(8+''-) + e-^(P'^+'''^))) (A.5) 



Looking at the second V'X diagram we have: 



c,i3,C 8q2(p-q) 



X [dABcdBCDie'^'''-^^-''^'''^ - e"^^ - 1 + e*(f-'^)'^) 



Now looking at ^XabXgh- 



c,d,e,/,B,c 32q2(p — q) 



(A.6) 



X I (e 

X 
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That is 



J ^^""i <i,e,B,cl6q2(p-q) 

-{h^g) 

^ ' (i,_B, c 8q2(p - q) 

-fABcfBCD{e~'^^'~'^^'+''^'^ + 1 + e"*'^'' + e^(*'s-«(9+'i+'>)))) 

-{h^g) (A.7) 

(Note that we also need to take into account the diagram where the internal ■^x is flipped. 
It is the same as what we have but with a ■H- g, b h and p -H- —p. We may for 
convenience take q o — q. We pick up an additional minus sign in the fABcfsCD term 
due to the differing order of the group factors.) 



B. Vanishing of One Loop Fermion Propagator at Zero Momentum 

Starting with the first diagram and using the simplified forms of the integrals derived in 
Appendix A (assuming an IR regulator) we have: 



-fr?Vd(0) = 



E 

BC 



1 



8q2q2 



(1-e- 



2^fABcfBCD{l + e"''^) 



= 



(B.l) 



as (1 — e '^<^)(1 + e'^<*) = 2i sin q^i and then the integrand is the combination of an odd and 
an even function. Next we calculate: 



E E s^2?;2 



d\ 



-iqm 



1) 



X {SagSmh + Sah^md^dg + ^ad^mg^dh — ^ah^rag " ^ag^md^dh — ^ad^mh^dg) 

x{dABcdBCD{e''^^'+''^ - 1 - e^'?(9+''-'^) + e"^^") 
-fABcfBCDie"''-^+'^^ + 1 + e^^(^+''-) + e-^^")) 



/ 



^ ' rn B.C 



x{dABcdBCD{e"^'^''+'''^ - 1 - e"^(3+h-a) + ^-iqa^ 
-fABcfBCD{e"''-o+'^^ + 1 + e^5(^+'>-) + e-''^-)) 



(B.2) 
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Then we can use the fact that if a ^ g and a ^ h then the expression disappears. If 
a = g = h again the expression disappears. So assuming a = h and a j/^ g we get: 



-fABcfBCD{e"'^^+'^^ + 1 + e"^^ + e-"^'^)) 

^^2^2 [dABcdBCD (sin + sin qg - sin ) 

-fABcfBCoismqa - sinQg - sing(„+g))) 
= (B.3) 



which vanishes term by term. We then move onto the second ipx diagram. 



c,B,C 

X ( dABcdBCD{e'''(^+''^ - - 1 + e-^^") 



-fABcfBCD{e''i^'^+'^) + e^«= + 1 + e-^««) (B.4) 



In a similar way to the previous diagram if a 7^ d and a 7^ e then the diagram vanishes. If 
a = d = e it also vanishes so we only need to deal with the case a = d, a ^ e: 



■^qa 



B,C 

-fABcfBCD{e"'^^+'^^ + e"i^ + 1 + e-'^- 

J X! g^2^2 (^""'^ ~ ^) (dABcdBCD{sinqa + sing(„+e) + singe) 

-fABcfBCoisinqa + sing(„+e) - sin^e) 
(B.5) 
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which again vanishes term by term. Finahy we show that Ixxi^) ~ ^■ 




X {dABcdBCD{e''^^'+''^ - 1 - e~^«<^ + e-^«(9+<^+'') ) 
-fABcfBCD{e''i^<'+'') + 1 + e-^^rf + e-^«(s+d+''))) 

f d'^q ^ 2i 



X (dABcrffiCD (sin - sing^ - sing(g+;j)) 

-fABcfBCD{smq(^d+g+h) + sm - sing(g+^))) 
= . 



(B.6) 



C. Coupling Constant Independence in A/^ = 4 SYM 

The continuum twisted J\f = 4 SYM discussed in this paper possesses a privileged set of 
operators whose expectation values can be shown to be independent of the background 
metric and hence topological. The condition for this to be true is that the operator be 
annihilated by the charge Q. In addition the expectation values of these operators can be 
shown to be independent of the coupling constant. As we will see this property remains 
true in the lattice theory and provides powerful constraints on the renormalization of such 
operators. To see this result consider the twisted lattice action which is the sum of Q- 
exact and Q-closed terms. The coupling constant dependence of the Q-closed term can 
be removed, without disturbing the Q BRST transformation, by rescaling the fields in 
appropriate ways. We show this below. 
The twisted action is 



9 ^exact ~1~ 9 ^closed 

9 9 





(C.l) 



A simple rescaling of the fields 



Xab Xab/9, #a, V V/9 



(C.2) 



gives the action 





9- 



(C.3) 
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Calling /3 = and writing the action as S = QA + iSdosed the expression for the 
expectation value of a Q-invariant operator O becomes 

Differentiating this expression with respect to /3 leads to 

= {Qlv)p{0)p-{OQK)p 
= {QA)p{O)0 - {Q{OK))p 

= , (C.5) 

where we have used the fact that as long as the BRST symmetry is not broken sponta- 
neously, the expectation value of the Q variation of some operator vanishes. Thus expec- 
tation values of Q-invariant observables are independent of (3 and hence can be computed 
exactly in the semi-classical limit (3 ^ oo. In this limit we need only do one loop calcula- 
tions around the classical vacua. 
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